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Abstract 

Conventional cone jet algorithms arose from heuristic considerations of LO hard scattering cou- 
pled to independent showering. These algorithms implicitly assume that the final states of individ- 
ual events can be mapped onto a unique set of jets that are in turn associated with a unique set of 
underlying hard scattering partons. Thus each final state hadron is assigned to a unique underly- 
ing parton. The Jet Energy Flow (JEF) analysis described here does not make such assumptions. 
The final states of individual events are instead described in terms of flow distributions of hadronic 
energy. Quantities of physical interest are constructed from the energy flow distribution summed 
over all events. The resulting analysis is less sensitive to higher order perturbative corrections and 
the impact of showering and hadronization than the standard cone algorithms. 
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A. Introduction 

The goal of a jet algorithm is to provide a precise mapping between the observed, long 
distance hadronic final states in high energy interactions and the underlying energetic par- 
tons participating in the true short-distance, hard scattering process]]]]. To appreciate this 
goal imagine comparing the observed final states detected by "real" detectors that have 
sizes of order centimeters to meters to what would be observed with a detector whose size 
is characterized by a distance scale of a fraction of a fermi. Between these two scales, 
i.e., from a fraction of a fermi to centimeters, the short-distance state evolves via higher 
order perturbative processes and the physics associated with showering and hadronization. 
These processes allow the short-distance partons, along with the spectators, to evolve into 
the observed hadrons. During the evolution the 4-momentum associated with an initial 
short-distance parton is spread out over a number of final state hadrons occupying an ex- 
tended region of phase space. To achieve our stated goal, the jet algorithm tries to identify 
these "related" hadrons into a single jet, whose total 4-momentum should track that of 
the initial parton. To perform this task with precision it is important that the results of 
applying the jet algorithm are insensitive to both higher order corrections and fluctuations 
in the showering/hadronization processes. The results of the jet algorithm should also be 
insensitive to the smearing effects of the detection process itself, e.g., due to the granularity 
of the detector. 

Current jet algorithms attempt to achieve the stated goal in a quite singular way by 
assigning the observed hadrons to unique jets on an event-by-event basis. This identification 
proceeds in the face of the fundamental fact that such a unique assignment cannot be precise. 
While the underlying parton that initiates the jet is a nonsinglet under the color symmetry 
of QCD, the hadrons are all color singlets. At the very least, the jets in the final state 
must represent the correlated evolution of more than one short-distance parton or of a 
short- distance parton correlated with a spectator parton, i.e., a color singlet set of initial 
partons. In some sense the most extreme approach is represented by cone algorithms as 
advocated at the 1990 Snowmass Workshop [0]. Cone algorithms associate hadrons into 
jets by identifying those that are nearby in angle. The underlying assumption is that the 
extra radiation produced by higher order corrections, showering and hadronization around 
an energetic parton appears symmetrically, i.e., occurs independently of the other color 
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charged objects in the final state. The alternative jet algorithm, the Kt algorithm ||, uses 
nearness in momentum space to identify the members of a jet and, at least to some extent, 
recognizes the "color-connectedness" of the radiation producing the final state. However, 
as noted above, both of these algorithms associate individual hadrons with unique jets on 
an event-by-event basis. We know that this procedure is only an approximation (due to 
the color conservation issue) and can lead to undesirable dependence (at least at the 10% 
level) on the details of the showering and hadronization processes. In the context of cone 
jet algorithms this latter point is discussed in more detail in another contribution to these 
Proceedings^. 

The Jet Energy Flow (JEF) approach described in this note is a simplified version of 
the more completely developed C-continuous observables or C-algebra formalism of F. V. 
TkachovQ for describing energy flow in hadronic collisions||. JEF accepts the reality 
that the hadronic final state represents the collective radiation from several out-flowing 
color charges, i.e., the underlying short- distance partons. No attempt is made to associate 
individual hadrons with unique jets, i.e., with unique underlying partons, on an event- 
by-event basis. Yet the energy flow pattern of an event still provides a footprint of the 
underlying partons, from which much of the same information provided by the standard 
jet algorithms can be extracted. As the subsequent discussion will indicate, it is a more 
reliable characterization of the event in the sense of exhibiting a reduced sensitivity to the 
showering and hadronization processes. The challenge in the JEF type analysis is to define 
observables that offer an informative comparison between theory and experiment. 

B. The JEF Formalism 

A primary strength of the JEF approach is that, in contrast with the usual algorithmic 
approach to jet identification, the JEF formalism generates, event-by-event, a smooth dis- 
tribution to characterize each event. In that sense, the JEF formalism is more analytic. For 
example, in the application of the cone algorithm the goal of identifying unique jets leads to 
the "stability" constraint A set of hadrons or partons that lie within a cone of a defined 
size R are identified as constituting a jet if and only if the energy-weighted centroid of the 
set of particles coincides with the geometric center of the cone. This constraint results 
in the non-analytic structure of the implementation of the algorithm, typically in the form 
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of step functions with complicated arguments. Only limited (and typically complicated) 
regions of the multi-particle phase space contribute to a jet. No such constraint arises in 
the JEF analysis. This distinction has several important consequences. 



1. The more inclusive and analytic calculations characteristic of a JEF analysis are more 
amenable to resummation techniques and power corrections analysis in perturbative 
calculations. 

2. Since the required multi-particle phase space integrations are largely unconstrained, 
i.e., more analytic, they are easier (and faster) to implement. Programs like JETRAD 
spend considerable computer time simulating the complicated phase space required by 
the algorithmic style jet algorithms. 

3. Since the analysis does not identify jets event-by-event, the analysis of the experimental 
data from an individual event should proceed more quickly. 

4. Signal to background optimization can now include the JEF parameters (and distri- 
butions). One cannot typically optimize a standard jet algorithm except for a limited 
number of parameters. 

We can define the fundamental distribution of the JEF analysis as follows. We start with 
a set of 4- vectors, p^ = (E,p x ,p y ,p z ), that represent either the partons in a perturbative 
calculation or the hadrons in a simulated or real event. In the latter case these 4- vectors 
might correspond as well to the location and energy deposited in individual calorimeter 
cells. If a given event corresponds to the measurement of N such 4- vectors, {p^l^Li = 
| (e\ P 1 ^ j , we have the 4-vector distribution for that event defined by 

N 



i=i 



p 



with the 2-dimensional 



where the directional unit vector is defined by P (m) — P / 
angular variable defined as m = (9, 0) (typical of lepton colliders) or m = (rj, 0) (typical 
of hadron colliders, where i] is the pseudorapidity, i] = In (cot 9/2)). This expression 
defines the underlying energy flow via E (m) = P (m) with a corresponding expression 
for the underlying longitudinal momentum flow P z (m). For the case of hadronic colliders 
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FIG. 1: Transverse energy (in GeV color coding) calorimeter readout in a typical CDF event as a 
function of pseudorapidity and azimuthal angle. The reconstructed CDF cone jets are indicated 
by the circles (with the same GeV color coding). 



the more familiar underlying transverse energy flow is defined by the composite quantity 
Et (m) = \jPx ( m ) + Py ( m ) or Pt ( m ) — E (m) x sin^ (m). Clearly many quantities can 
be constructed from the 4-momentum distribution of Eq. |TJ, including the usual cone jet 
algorithm. The Et (r], 0) distribution for a typical CDF jet event is illustrated in Figure [l] 
along with the cone jets "found" with the CDF cone jet algorithm. 

Using the underlying 4-vector distribution we define the jet energy flow (JEF) via a 
smearing or averaging function A as 

J M (m) = J dm' P M (mf) x A(m! - m) , (2) 

where A is normalized as 

J dm A (m) = 1. (3) 

A simple (but not unique) form for the averaging function in terms of the general 2-tuple 
of angular variables m = (a, (3), which provides a direct comparison with the jet cone 
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FIG. 2: The Ej> (77, 0) flow (including the factor irR 2 ) using the CDF event of Figure 1. 



algorithm, is 



A (m) = A (a, (3) 



Q(R-r(a,/3)) 



6 (r - ^/a 2 + (3 2 



(4) 



irR 2 ttR 2 

where R is the cone size and r (a, f3) is the distance measure in the space defined by (a, 0). 
For comparison with the existing jet cone analyses we will discuss the case m = (rj, <fi). As 
a specific example we exhibit the jet transverse energy flow (transverse JEF) 



m) 



(5) 



times a factor of ttR 2 (Et = ttR 2 x Jt) for the case of R = 0.7 in Figure and |3| for the 
same event displayed in Figure |I]. Clearly the same general structure is present in all three 
figures. Note that the transverse JEF is smeared on a scale R compared to the underlying 
Et distribution of Figure 0. For comparison the Snowmass cone jet algorithm^, identifies 
jets at a discrete set of values of locations rrij defined by the Et weighted cone "stability" 
constraint. These stable cone locations rrij are the solutions of the equation 



dm' Et (m 1 ) x (m! — rrij) x A [m — rrij) = 0. 



(6) 
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FIG. 3: Same as Figure 2 except that the maximum Et for the color coding is 15 GeV. 

The corresponding cone jet Et values are found by evaluating Eq. [5] (times irR 2 ) at the 
jet positions, Ej>j = nR 2 X Jt {trij). The non-analytic character of the jet cone algorithm 
referred to earlier arises from the need to solve Eq. ^ and then work with only the discrete 
set of solutions, i.e., the jets in an event. 

C. Observables 

We can now proceed to define more general observables. The basic assumption of the JEF 
approach is that event-by-event each value of the direction m is equally likely to correspond 
to a jet with 4-momentum proportional to J^{m). Relative probabilities of observables 
having values in a specified range will correspond to the size of the area in m covered by the 
JEF with the correct range of values. To illustrate these ideas, consider a general nth order 
observable C n represented by an nth order function of J M (m), C (J (mi) , . . . , J (m n )). The 
corresponding event probability distribution, including the possibility of providing a set of 
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angular cuts f2, is given by 

P(C n \n(m cut ))= (jlj ^n( mi - mcut )^6(C n -C(J( mi ),... ,J(m n ))) (7) 

(icTjEF 

dC n 

where we have normalized the area to the "cone size" nR 2 . (Note that one might also 
consider applying a cut directly on the J(m), e.g., Jo{m) > E cut .) To determine the 
differential cross section for the observable C n from an experiment we simply sum over 
events as 

C ^= E^i^k^)), ( 8 ) 

n events 

where £ is the integrated luminosity. We obtain an event occupancy probability O by 
integrating over the probability function 

»C„(max) 
' C n (min) 



0(C n (min),C n (max)|fi(m cut )) = / dC n P (C n \Q (m cut ) ) (9) 

J C n (min) 

= j ~^2^ ^ ~~ mcut ^j ( Cn ( maX ) ~ C ( J ( m l) > • • ■ > J i m n))) 

x (C ( J (mi) , . . . , J (m n )) - C n (min)) . 

This final expression indicates that we are simply calculating the relative area in m for which 
the JEF has the correct value to yield the desired value of the observable. We can then 
count the number of events with effective occupancy number O and convert it into a cross 
section, 

Ca (C n (min) , C n (max)) = ^ O (C n (min) , C n (max)| Q (m cut )) . (10) 

events 

This formula can be used to construct bin values and the corresponding distribution. Let 
us illustrate these ideas by considering some explicit examples. 



1. The JEF jet mass M ( J (m)) = irR 2 x a/ J u (m) J> (m) is an example of a C\ observable 
with a event probability distribution of the form 

P(Mj) = j ^ HMj - M (J (m))) . (11) 
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FIG. 4: The JEF jet invariant mass as denned in Eq. 11 for W boson decay in Higgs boson decay. 



The corresponding occupancy probability has the form 

/(JTff 
9 (Mj (max) — M (J (m))) 9 (M (J (m)) - Mj (min)) 

Thus the fraction of the events with a JEF jet with mass in the specified range is pro- 
portional to the fractional area in the m plane occupied by JEF jets with a mass value 
in that range. To obtain the final distribution we sum over events. The (simulated) 
JEF jet mass distribution for the W decay into hadrons (treated as a single jet) in the 
process pp — > H + X — > W + W~ + X — > Iv + hadrons + X is exhibited in Figure [|. 

2. The JEF jet transverse energy Et in the variables appropriate to a hadron collider is 
another C% observable. The relative probability distribution for a CDF type rapidity 
acceptance and CDF Et definition looks like 

P (E T \Q (0.1 < \ V \ < 0.7) ) = y°' 7 d \v\ j d<p5 (E T -E{J ( V , 0)) x sin (9 (??))) , 

(12) 

where the JEF energy distribution is given by E (J (77, 0)) = nR 2 x Jo(i],<j)). As 
suggested above, we can obtain the effective occupancy number of JEF jets (per event) 
above an energy cut E T , m m by integrating 

O(E Ttmin \tt(0.1 < \rj\ < 0.7)) = / dE T P(E T \tt(0.1 < \r]\ < 0.7)). (13) 



9 




20.0 40.0 60.0 80.0 100.0 




1 ' T~ 

).0 20.0 40.0 60.0 80.0 100.0 



4.0 — 

3.5 - 
3.0 - 
2.5 - 
2.0 - 
1.5 - 
1.0 - 
0.5 - 
0.0 — 



0.0 20.0 40.0 60.0 80.0 100.0 



FIG. 5: The Et probability function of Eq. 12 for the event of Figure 1: the left figure is for a bin 
width of 1 GeV and middle figure is for a bin width of 5 GeV. The right figure is the occupancy 



number of Eq. 13. The individual data points corresponding to the 5 largest energy jets found by 



the CDF cone jet algorithm and illustrated in Figure 1. 



These quantities as evaluated for the sample jet event of Figure 1 are illustrated in 
Figure [|. The jets found by the standard CDF cone jet algorithm are also indicated 
as data points in the figures and correlate well with the peaks in the JEF probability 
distribution. 

3. The JEF di-jet invariant mass M (J, J) is an example of a C 2 observable with the 
form M 2 (J (mi) , J (m 2 )) = (vri? 2 ) 2 x (J(m 1 ) + J(m 2 )) fi (J(m 1 ) + J(m 2 )) 11 , which 
assumes that the two JEF jets are non-overlapping.. The corresponding probability 
distribution is defined by 

P (M 2 jj) =U (M'j - M 2 (J ( mi ) , J (ma))) . (14) 

D. An Example JEF Analysis 

As an example of a JEF style analysis we briefly review a JEF di-jet analysis performed 
previously ||. The goal is to calculate the differential transverse energy distribution of a 
jet in the CDF central rapidity strip, 0.1 < \r]\ < 0.7, while requiring that a second jet with 
transverse energy at least as large as one half of the transverse energy of the central jet, 
Et,2 > -Et,i/2, is tagged in a forward region, 1.2 < j^l < 1.6. The corresponding di-JEF 
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FIG. 6: Differential cross section of Eq. ^ for both JEF (moving cone) and EKS (fixed cone) 
analysis at LO and NLO. 



Pdi-jct (E T ) = —^—2 1 1 1 1 d^dfad^dfc n (0.1 < \rjl\ < 0.7) x n (1.2 < M < 1-6) 




probability density function for an event obeying the appropriate cuts is expressed as 
1 

(15) 

x Q (E T (J ( % , 2 )) > E T (J ( Vl , 01)) /2) x5(E T - E T (J ( Vl , 00)) . 
Using this probability function the desired differential cross section is obtained from Eq. H 

events 

The perturbative results for this cross section at LO and NLO for both the JEF analysis 
of Eq. [16] and using the standard cone jet analysis of EKS0] are exhibited in Figure 6 
for the case R = 0.7. These results illustrate some of the desirable features of the JEF 
approach. Note first that at smaller E? values both jet definitions yield similar results at 
both LO and NLO. However, at larger values, where the boundaries of phase space 
play are more relevant, the JEF result is larger than the "traditional" cone jet result and, 
more importantly, is less sensitive to the higher order corrections. We can understand this 
reduced sensitivity in terms of the smearing of the rapidity cuts in the JEF analysis. In the 
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JEF analysis the underlying partons can violate the rapidity cuts by as much as R and still 
contribute to a JEF style jet that respects the cut. For example, the parton contributing 
to the secondary jet is only required to obey 1 772 1 > 1.2 — R = 0.5 in order to make a nonzero 
contribution. In contrast, the traditional cone jet with a single parton inside requires the 
parton to be collinear with the jet and thus in the current analysis receives contributions 
only from j^l > 1-2. This smearing of the details of the rapidity cuts explains both the 
larger magnitude and the reduced dependence on higher orders of the JEF di-jet analysis. 
We can expect a similarly reduced dependence on the stochastic effects of showering and 
hadronization. 

E. Concluding Remarks 

We have discussed a different approach to the jet analysis of hadronic states, the JEF 
analysis, which follows from the earlier C-algebra formalism 0] and which differs from tradi- 
tional algorithmic approaches in that unique jets are not identified event-by-event. Instead 
the analysis proceeds through the evaluation of Jet Energy Flow distributions. The brief 
discussion presented here suggests that the JEF style analysis of hadronic final states in hard 
scattering processes will provide observable measures of the underlying short-distance parton 
structure that are less sensitive to higher order corrections and to showering/hadronization 
corrections than more conventional jet algorithm analyses. Clearly much more needs to be 
done in order to demonstrate and make use of this conclusion. 
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